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Abstract 

We revisit the vertical structure of black hole accretion disks in spherical coordi- 
nates. By comparing the advective cooling with the viscous heating, we show that 
advection-dominated disks are geometrically thick, i.e., with the half-opening angle 
A6' > 27r/5, rather than slim as supposed previously in the literature. 
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1. Introduction 

It was known long since that the very basic assumption of the Shakura-Sunyaev disk 
(SSD, Shakura & Sunyaev 1973), that is, the geometrical thinness of the disk, H/R<^ 1, where 
H is the half thickness of the disk and R is the radius in cylindrical coordinates, would break 
down for the inner region of the disk in some specific situations. For example, when the mass 
accretion rate M approaches and surpasses its critical value corresponding to the Eddington 
luminosity, radiation pressure will act to huff the inner region of the disk in the vertical direction; 
or when the cooling mechanism is inefficient, so that the temperature in the disk becomes very 
high, then gas pressure will act in a similar way. In either of these two situations, the inner 
region of the disk will get geometrically thick, i.e., with H/R ~ 1 (e.g., Frank et al. 2002, 
p. 98). Based on these understandings, two types of models were proposed more than twenty 
years ago, namely the optically thick, radiation pressure-supported thick disk (Abramowicz et 
al. 1978; Paczyhski & Wiita 1980; Madau 1988) and the optically thin, ion pressure-supported 
thick disk (Rees et al. 1982). To avoid mathematical difficulties, in these models the disk was 
assumed to be purely rotating, i.e., with no mass accretion. However, the very existence of non- 
accreting thick disks was thrown into doubt by the discovery of Papaloizou & Pringle (1984) 
that such disks are dynamically unstable to global non-axisymmetric modes. Since the work of 
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Blaes (1987), it had been recognized that it is accretion, i.e., radial matter motion and energy 
advection into the central black hole, that can sufficiently stabilize all modes. Accordingly, the 
concept of advection dominance was introduced and two new types of models were constructed, 
namely the optically thick, radiation pressure-supported slim disk (Abramowicz et al. 1988) 
and the optically thin, ion pressure-supported, advection-dominated accretion flow (ADAF, 
Narayan & Yi 1994; Abramowicz et al. 1995). Both these two types of models are popular 
nowadays. 

Slim disks and ADAFs were supposed to be geometrically slim, i.e., with H/W^l, neither 
thin nor thick. The reason for this restriction is the following. As argued by Abramowicz et 
al. (1995), the advection factor /adv = <5adv/Qvis, where Qadv is the advective cooling rate per 
unit area and Qvis is the viscous heating rate per unit area, should satisfy the relation 

/adv>ffV. (1) 



Obviously, advection can be important only for disks that are not thin. But the disk cannot 
be thick either, because the value of /adv cannot exceed 1. 

Recently, Gu & Lu (2007, hereafter GL07) addressed a problem in the slim disk model of 
Abramowicz et al. (1988, see also Kato et al. 1998). In this model, the gravitational potential 
was approximated in the form suggested by Hoshi (1977), i.e., 

iPiR,z)c^ij{R,0) + ^nlz^ , (2) 

where fix is the Keplerian angular velocity. As shown by GL07, such an approximation is valid 
only for geometrically thin disks with H/R ^ 0.2, and for a larger thickness it would greatly 
magnify the gravitational force in the vertical direction. Accordingly, the widely adopted 
relationship HQi^/cs = constant can approximately hold only for thin disks as well. Since 
formula (1) was derived by using this relationship, its validity for thicker disks has not been 
justifled. GL07 noted that, when the vertical gravitational force is correctly calculated with 
the explicit potential ip{R,z), "slim" disks are much thicker than previously thought. However, 
the work of GL07 was still within the framework of the slim disk model in some sense. In 
particular, those authors did not consider the vertical distribution of velocities, but instead 
kept the assumption of vertical hydrostatic equilibrium, 

1 dp dip ^ 

pdz dz 

which is a simpliflcation of the more general vertical momentum equation 

Idp ^ dip ^ dvz ^ dvz Q 
pdz dz ^ dR ^ dz 

(e.g., Abramowicz et al. 1997), where p is the mass density, p is the pressure, and vr and 
Vz are the cylindrical radial and vertical velocities, respectively. While the terms containing 
Vz in equation (4) can be reasonably dropped for thin disks because in this case Vz must be 
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negligibly small, it needs a careful consideration whether the same can be done for not thin 
disks (Abramowicz et al. 1997, also see below in §2). 

Also regarding to the two main features of advection-dominated disks, i.e., the advection 
dominance and the slimness, an important different approach was made earlier by Narayan & 
Yi (1995, hereafter NY95). NY95 considered rotating spherical accretion flows ranging from 
the equatorial plane to the rotation axis, i.e., with H/ R oo and with no free surfaces. They 
assumed self- similarity in the radial direction and solved differential equations describing the 
vertical structure of the flow, and showed that, comparing to their exact solutions, the solu- 
tions obtained previously with the vertical integration approach are very good approximations, 
provided "vertical" means the spherical polar angle 9, rather than the cylindrical height z. 
This seemed to indicate that advection-dominated disks are not necessarily limited to be slim. 
However, those authors did not calculate the advection factor /^^^ (they defined fl^^ = q'adv/Q'vis, 
with gadv and gvis being the advective cooling rate and the viscous heating rate per unit volume, 
respectively), but rather set it a priori to be a constant. It is still not answered how their f'^^ 
varies with 6', or how /adv per unit area varies with the thickness of the disk, and what is 
required for advection to be dominant. 

In this work we try to make some complementarity to NY95 and some refinements to 
GL07. We consider the vertical structure of accretion flows with free surfaces and show that 
advection-dominated disks must be geometrically thick rather than slim. Our results may 
suggest to recall the historical thick disk models mentioned above, but with improvements that 
they have to include accretion now. 

2. Equations 

We consider a steady state axisymmetric accretion flow in spherical coordinates (r, 
9, 0) and use the Newtonian potential if) = —GM/r since it is convenient for the self-similar 
formalization adopted below, where M is the black hole mass. The basic equations of continuity 
and momenta are 

Id Id 

— — (rVr) + —^JIB^^^'^^P'^^o) = , (5) 
or r sm 9 o9 

^ ^ ^ , ^ ^'j. ^ , GM 1 dp . , 

^ dr r \d9 I r pdr ^ 




^ dr r \d9 ^ ) r pr d9 ^ 

dr r d9 r pr-^ dr 

(e.g., Xue & Wang 2005), where Vr, vg, and are the three velocity components. We 
assume that only the r0-component of the viscous stress tensor is important, which is 
tr(p = i^prd{vtp/r)/dr, where u = acj.r/vK is the kinematic viscosity coefficient, a is the con- 
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stant viscosity parameter, Cs is the sound speed defined as cj. =p/ p, and vk = [GMjrY^'^ is the 
Keplerian velocity. 

We do not simply assume vertical hydrostatic equilibrium (eq. [3]). Equation (7) is the 
general vertical momentum equation in spherical coordinates, corresponding to equation (4) 
in cylindrical coordinates. Abramowicz et al. (1997) have given several reasons why spherical 
coordinates are a much better choice. We only mention one of these reasons that is partic- 
ularly important for our study here. The stationary accretion disks calculated in realistic 
two-dimensional (2D) and three-dimensional (3D) simulations resemble quasi-spherical fiows, 
i.e., in spherical coordinates the half-opening angle of the fiow A6' ^ constant, or in cylindrical 
coordinates the relative thickness H/ Rk. constant, much more than quasi-horizontal fiows, i.e., 
H K, constant (e.g., Papaloizou & Szuszkiewicz 1994; NY95). If no outfiow production from the 
surface of the disk is assumed, then obviously f ^ = is a reasonable approximation for disks 
with any thickness (Xue & Wang 2005); but cannot be neglected for not thin disks because 
there is a relation Vz/vr ~ H/ R for quasi-spherical fiows, making equation (4) difficult to deal 
with. 

Similar to NY95, we assume self-similarity in the radial direction 

Vr (X r~"^/^; ve = 0; v^p oc r'^^"^; 

p cx r~'^/^; Cs oc r~^^'^. 

The above relation automatically satisfies the continuity equation (5). By substituting the 
relation, the momentum equations (6-8) are reduced to be 

i«? + fc,^ + «J-4 = 0, (9) 
fi-4co.O, (10) 

Four unknown quantities, namely Vr, v^, Cg and p, appear in these three equations. This is 
because we do not write the energy equation, whose general form is gvis = 9adv + Imd, where 
grad is the radiative cooling rate per unit volume. In principle, the general energy equation 
should be solved, and then /^^^ is obtained as a variable, as done, e.g., by Manmoto et al. 
(1997) for ADAFs and by Abramowicz et al. (1988) and Watarai et al. (2000) for slim disks. 
But due to complications in calculating the radiation processes, in NY95 and even in works 
on global ADAF solutions (e.g., Narayan et al. 1997), gadv = fUvQ^is or Qadv = /advQvis was 
used instead as an energy equation and f!^^^ or /adv was given as a constant. Since our purpose 
here is to investigate the variation of /adv with the thickness of the disk, we wish to calculate 
(5 adv and Qvis respectively, and then estimate /adv To do this, we further assume a polytropic 
relation, p = Kp'^, in the vertical direction, which is often adopted in the vertically integrated 
models of geometrically slim disks (e.g., Kato et al. 1998, p. 241). We admit that the poljd;ropic 
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assumption is a simple way to close the system, and then enables us to calculate the dynamical 
quantities and evaluate /adv self-consistently. 

With the polytropic relation and the definition of the sound speed = p/ p, equation 
(10) becomes 

_^.I_„JeotO, (12) 

which along with equations (9) and (11) can be solved for Vr, v^, and c^. A boundary condition 

is required for solving the differential equation (12), which is set to be = (accordingly p = 

and p = 0) at the surface of the disk. The quantities gadv = pVr{dlnp/dr — '-/din p/dr)/{'~f — 1) 

and gvis = i'pr'^[d{v^/r)/dr]'^ are expressed in the self-similar formalism as 

5-37 pvr , , 

5adv--^(^— , (13) 

(14) 

then Qadv and Qvis are given by the vertical integration, 
/-f+Ae 

Qadv = / gadv rsinO dO , (15) 
/•f+Ae 

Qvis = / gvis r^AnO de , (16) 

and /adv = Qadv/Qvis IS obtained. In our calculations a = 0.1 is fixed. 

3. Numerical results 

We first study the variation of dynamical quantities with the polar angle 6 for a given 
disk's half-opening angle A^. Figure 1 shows the profiles of Vr (the dashed fine), (the dot- 
dashed line), Cs (the solid line), and p (the dotted line) for three pairs of parameters, i.e., 7 = 4/3 
and = 0.257r for Fig. la, 7 = 4/3 and A^ = 0.457r for Fig. 16, and 7 = 1.65 and A^ = 0.4987r 
for Fig. Ic. The parameters are marked in Figure 3 by filled stars, which clearly show the 
corresponding values of the advection factor /adv Obviously, advection is not significant for 
case a (/adv < 0.1), but is dominant for cases b and c (0.5 < /adv < 1)- Comparing our results 
with Fig. 1 of NY95, it is seen that the profiles of Vj. and p are similar, i.e., Vr (the absolute 
value) and p increase with increasing 6 and achieve the maximal value at the equatorial plane 
{6 = vr/2). On the contrary, the two profiles of Cg are significantly different. In their Fig. 1, the 
value of Cs decreases with increasing 6 and achieves the minimal value at the equatorial plane; 
in our Fig. 1, however, Cg increases with increasing 6 and achieves the maximal value at the 
equatorial plane. In our opinion, the difference results from different assumptions, i.e., NY95 
assumed an energy advection factor /^^^ in advance, whereas we solve for the energy advection 
factor /adv self-consistently based on a polytropic relation in the vertical direction. We think 
that our profile for is reasonable for disk-like accretion. For example, in the standard thin 
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disk, the direction of the radiative flux is from the equatorial plane to the surface, which means 
that the temperature (or the sound speed) decreases from the equatorial plane to the surface. 
Such a picture agrees with our Fig. 1 but conflicts with Fig. 1 of NY95. 

Figure 2 shows the variation of /adv with for the ratio of specific heats 7 = 4/3. 
Advection dominance means 0.5 < /adv < 1- We first explain the two dashed lines and the 
dotted line that correspond to previous works in the slim disk model, then the solid line that 
represents our results here, and leave the dot-dashed line later. Both the two dashed lines are 
obtained by assuming vertical hydrostatic equilibrium (eq. [3]) and using the Hoshi form of 
potential (eq. [2]), thus the relation HQi^/cs = constant is adopted. The difference between 
these two lines is the following. For line a, the simple one-zone treatment in the vertical 
direction is made as in the SSD model; then in equation (3), dp/dz ~ —p/H, dip/dz ~ ^\^H^ 
and HQi^/cs = 1 is obtained (e.g., Kato et al. 1998, p. 80). For line b, there is some improvement 
in the sense that the vertical structure of the disk is considered. By assuming a polytropic 
relation, the vertical integration of equation (3) gives HQyi/cs = 3 (e.g., Kato et al. 1998, 
p. 242). Because of these different treatments in the vertical direction, these two lines show 
different variations of /adv with A6 and different maximum values of A6. The upper limit 
of /adv is 1 (full advection dominance), beyond which there would be no thermal equilibrium 
solutions. It can be analytically derived that for the case of line a, the maximum value of 
corresponding to /adv = 1 is A6'niax = arctan(y^2/7), or in cylindrical coordinates the maximum 
relative thickness (H/R)^^^ = y^2/7; and for the case of line b it is A6'max = arctan(3/2) or 
(if/i?)max = 3/2. As mentioned in §1, the thickness of the disk in the slim disk model had 
been underestimated because the vertical gravitational force was overestimated by the Hoshi 
form of potential. Even so, according to the more sophisticated version of the slim disk model 
(line b), advection dominance /adv > 0.5 would require H/R> 1 {A9 > 7r/4), and full advection 
dominance would require H/R = 3/2, in contradiction with H/R ^ 1, the supposed feature of 
the model. 

The dotted line in Figure 2 is for the results of GL07. The point made in that work 
was that the explicit potential iIj{R,z), rather than its Hoshi approximation (eq. [2]), was 
used, so that the vertical gravitational force was correctly calculated. But GL07 still kept the 
assumption of vertical hydrostatic equilibrium (eq. [3]), i.e., the terms containing in equation 
(4) were incorrectly ignored. Because of this, the thickness of the disk was overestimated; and 
accordingly, it seemed that advection dominance can never be possible, since even for the 
extreme thickness A^ = 7i/2 (or H/R 00) the value of /adv can only marginally reach to 0.5. 

We make improvements over GL07. We use spherical coordinates with the assumption 
vg = 0, which is better than = in cylindrical coordinates; and then calculate the vertical 
distribution of velocities {vr and Vcf,) and thermal quantities (p, p, and c^). Our results are 
shown by the solid line in Figure 2. It is seen that advection dominance (/adv > 0.5) is possible, 
but only for A6 > 27r/5 (or 72°). Therefore, advection-dominated disks must be geometrically 
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thick, rather than shm as previously supposed. 

It is also seen that line b, the dotted line, and the solid line in Figure 2 almost coincide 
with each other for thin disks with A6 ^ O.Itt. This is natural, since for thin disks both the 
Hoshi approximation of potential and the assumption of vertical hydrostatic equilibrium are 
valid, and the three approaches represented by the three lines make no significant difference. 
But the one-zone treatment, i.e., total ignorance of the vertical structure of the disk, seems to 
be too crude, making the resulting line a deviate from the other three lines even for thin disks. 

The value 7 = 4/3 in Figure 2 corresponds to the optically thick and radiation pressure- 
dominated case, to which the historical radiation pressure-supported thick disk and the slim 
disk belong; while it is 7 ^ 5/3 for the optically thin and gas pressure-dominated case, to which 
the historical ion pressure-supported thick disk and the ADAF belong. In Figure 3, the four 
solid lines show variations of A6 with 7 for four given values of /adv It is seen that advection 
dominance (/adv > 0.5) requires to be large for any value of 7; and that for a fixed /adv (the 
same degree of advection), the required AO increases with increasing 7, that is, for advection 
to be dominant, optically thin disks must get even geometrically thicker than optically thick 
ones. 

For the geometrically thin case, A6 <^ 1, the Taylor expansion of equations (9), (11), 
and (12) with respect to can be performed, and we derive an approximate analytic relation: 

^ (5-37)(27-l) 2 .^7^ 
37(57-3) ' 

which is similar to equation (1) in cylindrical coordinates. The dot-dashed lines in Figures 2 
and 3 correspond to equation (17) for a fixed 7 = 4/3 and for a fixed /adv = 0.01, respectively. 
It is seen from Figure 2 that, as expected, the analytic approximation of equation (17) agrees 
well with the correct numerical results (the solid line) for small A6, but deviates a lot for large 
A6. In Figure 3 a good agreement between equation (17) and the numerical results (the lowest 
solid line) is seen again, especially for small values of 7. The limitation that equation (17) is 
valid only for small A^, and accordingly only for small /adv, should also apply to equation (1), 
because that equation is derived with the Hoshi form of potential. 

4. Discussion 

The key concept of the slim and ADAF disk models is advection dominance. This concept 
was introduced rather as an assumption, whether and under what physical conditions can it be 
realized have not been clarified. The main result of our work is to have shown that, in order 
for advection to be dominant, the disk must be geometrically thick with the half-opening angle 
A^ > 27r/5, rather than slim as suggested previously in the slim disk and ADAF models. Thus, 
advection-dominated disks are geometrically similar to the historical thick disks metioned in 
§1. This result is obvious because, as revealed in GL07, in the slim disk and ADAF models the 
vertical gravitational force was overestimated by using the Hoshi's approximate potential, and 
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accordingly the disk's thickness was underestimated. NY95 considered accretion flows with no 
free surfaces and found that when the given advective factor /adv(=Q'adv/9vis) 1 (full advection 
dominance), their solutions approach nearly spherical accretion. If "nearly spherical" can be 
regarded as extremely thick, then their results and ours agree with each other, but we take 
a different approach. We do not give the value of /adv(= Qadv/Qvis) in advance, but instead 
consider accretion flows with free surfaces, i.e., accretion disks. The boundary condition is 
set to be p = 0, which is usually adopted in the literature (e.g., Kato et al. 1998). Then the 
thickness of the disk, A6, makes sense, and we calculate /adv to see how it relates to A6. 

Many 2D and 3D numerical simulations of viscous radiatively inefficient accretion flows 
(RIAFs) revealed the existence of convection-dominated accretion flows (CDAFs), while ADAFs 
could not be obtained (e.g.. Stone et al. 1999; Igumenshchev & Abramowicz 2000; McKinney & 
Gammie 2002; Igumenshchev et al. 2003). We think that this fact probably indicates that the 
existing analytic ADAF models might have hidden inconsistencies, and the incorrect treatment 
of the vertical structure might be one such inconsistency, as addressed in our work. Moreover, 
the recent radiation-MHD simulations (Ohsuga et al. 2009) showed that the disk is geomet- 
rically thick in their models A and C (corresponding to slim disks and ADAFs, respectively), 
which is in agreement with our results. 

Apart from the convective motion, the outflow is found in 2D and 3D MHD simulations of 
non-radiative accretion flows (e.g.. Stone & Pringle 2001; Hawley & Balbus 2002). For optically 
thick flows, the circular motion and the outflow are found in 2D radiation-HD simulations 
(e.g., Ohsuga et al. 2005; Ohsuga 2006). The assumption vg = would break down when 
the convective motion or the outflowing motion is significant, thus we have to point out the 
limitation of our solutions, which are based on the self-similar assumption in the radial direction 
and particularly for vg = 0. 

In this paper we have not shown the exact thermal equilibrium solution for a certain mass 
accretion rate. We wish to stress that our main concern here is the relationship between the 
energy advection factor and the thickness of the disk. The well-known formula (1), which was 
previously believed to be valid for both optically thick and thin disks, implied that advection- 
dominated accretion disks are geometrically slim. As shown in Figures 2 and 3, however, formula 
(1) is inaccurate for disks that are not geometrically thin. We think that the new relationship 
between /adv and A6, shown in Figures 2 and 3, should also work for both optically thick and 
thin cases. Even without the exact solutions, we can predict that advection-dominated accretion 
disks ought to be geometrically thick rather than slim. Our next work will concentrate on the 
optically thick disks and take the radiative cooling into consideration. In the vertical direction, 
we will solve the dynamical equations combined with the radiative transfer equations, thus the 
polytropic assumption will be relaxed. At that step, we will be able to calculate the thermal 
equilibrium solutions with given mass accretion rates and show the optical depth, pressure, and 
luminosity of the disks. 
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Jt/4 



71/2 



Fig. 1. Variations of Wr, v^, Cs, and p with the polar angle 9 for three pairs of parameters: (a) 7 = 4/3 
and A9 = 0.257r; (b) 7 = 4/3 and A6I = 0.457r; (c) 7 1.65 and AO = 0.4987r. 
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Fig. 2. Variation of the advection factor /adv with the disk's half-opening angle for the ratio of 
specific heats 7 = 4/3. The solid line shows our numerical results. The dot-dashed line corresponds to the 
analytic approximation of equation (17). The two dashed lines are for the previous results in the slim disk 
model with the Hoshi form of potential, and the dotted line is for the previous results in GL07. 



11 




Y 



Fig. 3. Variation of A6' with 7 for given values of /adv The sohd lines show numerical results, and 
the dot-dashed line corresponds to equation (17). The three filled stars denote the parameters chosen in 
Figure 1. 
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